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DERIVATION  OF  A  WAVENUMBER  FILTER  DESIGN  FOR  THE 
MEASUREMENT  OF  TURBULENT  WALL  PRESSURE  FLUCTUATIONS 


INTRODUCTION 


The  measurement  and  analysis  of  turbulent  wail  pressure  fluctuations  are  important  to  problems 
of  flow-induced  vibration  and  noise  as  well  as  to  studies  of  turbulence  physics.  The  space  time 
Field  p(xi,X2,t)  and  associated  three-dimensional  wavenumber-frequency  spectrum  d>(ki,k2,co)  are 
of  primary  interest.  To  date,  they  have  rarely  been  directly  measured  in  laboratories  due  to  the 
extensive  number  of  sensors  and  signal  processing  required.  A  common  approach  used  recently 
by  Karangelen  et  al.  (1991)  and  Keith  and  Barclay  (1993)  has  been  to  estimate  the  wavenumber- 
frequency  spectrum  from  cross-spectral  measurements  made  with  a  small  number  of  sensors  at 
unequal  spacings.  Although  this  approach  has  been  useful  for  estimating  the  convective  ridge 
region,  uncertainties  arise  due  to  the  limited  number  of  spatial  separations  at  which  data  arc 
obtained.  In  particular,  these  uncertainties  preclude  accurately  e  d mating  the  lower  spectral  levels 
at  wavenumbers  outside  of  the  convective  ridge  region. 

Many  single  point  experimental  measurements  of  p(t)  and  the  one-dimensional  spectrum  d>(co) 
have  appeared  in  the  literature,  as  reviewed  by  Keith  et  al.  (1992).  Here,  we  determine  a 
wavenumber  filter  design  for  the  measurement  of  the  two-dimensional  spectrum  C>(ki,to). 

Although  arrays  of  sensors  are  commonly  utilized  for  acoustic  measurements,  they  are  generally 
not  implemented  in  incompressible  fluid  mechanics  experiments.  The  methodology  is  also 
applicable  for  fluctuating  wall  shear  stress  measurements.  However,  a  widely  accepted  model  for 
the  wall  shear  stress  wavenumber-frequency  spectrum  required  in  this  analysis  does  not  exist  at  the 
present  time.  This  is  largely  due  to  the  lack  of  published  research  on  correlation  or  cross-spectral 
density  measurements  of  the  wall  shear  stress  field. 

We  are  primarily  concerned  with  three  array  parameters:  size,  spacing,  and  number  of  sensors. 
Following  Blake  and  Chase  (1971),  Stnol'y&kov  aruj  Tkachenko  (1983),  Farabee  and  Geib 
(1991),  and  Manoha  (1991),  we  use  a  linear  array  of  equally  spaced  sensors.  The  mechanical  and 
electrical  design  of  the  sensors  will  not  be  addressed  here.  Currently,  the  size  of  available  sensors 
poses  limitations  on  the  design  of  a  wavenumber  filter.  Here  we  assume  that  microfabrication 
technology  will  allow  a  nearly  optimum  design  to  be  specified  at  the  outset.  The  design  problem 
then  involves  estimating  the  true  field  to  be  measured  and  determining  the  effect  of  the  above  three 
parameters  on  the  measured  field.  The  sensors  are  assumed  to  have  a  constant  frequency  response 


and  uniform  sensitivity  over  the  sensing  area.  We  first  briefly  review-past  investigations  of  sensor 
size  effects  on  the  wall  pressure  frequency  spectrum  d>(«). 

THE  EFFECT  OF  SENSOR  SIZE  ON  THE 
WALL  PRESSURE  FREQUENCY  SPECTRUM 

The  early  work  of  Corcos  (1963)  showed  that  an  increase  in  sensor  size  leads  to  a  more  rapid 
roll-off  of  the  measured  frequency  spectrum  <J>(co)  at  high  frequencies.  This  effect  has  been  well 
established  and  is  supported  with  the  comparisons  of  published  measurements  of  <I>(W)  by  Keith  et 
al.  (1992).  Schewe  (1983)  measured  wall  pressure  spectra  with  circular  sensors  of  varying 
diameter  and  concluded  that  the  diameter  d+  =  dui/v  =  19  was  adequate  to  resolve  O(co)  at  the 
highest  energy -containing  frequencies  in  the  turbulent  field.  Schewe  did  not  report  measurements 
for  d+  <  19.  The  comparisons  presented  by  Keith  et  al.  (1992)  showed  Schewe's  spectral  levels 
(for  d+  =  19)  scaled  with  inner  variables  compared  well  with  the  numerical  simulation  of  Choi  and 
Moin  (1990)  (for  d+~  18).  However,  Keith  et  al.  (1992)  also  showed  that  the  Corcos  (1963) 
correction  increased  Schewe's  spectral  levels  (ford*  =  19)  by  2.5  dB  at  the  highest  frequencies 
measured. 

The  pressure  sensor  diameter  d+  required  for  proper  resolution  of  0(a>)  at  all  frequencies  and 
Reynolds  numbers  remains  unknown.  Schewe's  results  showed  that  sensors  with  d+  >  19  are 
inadequate.  Head  and  Bandyopadhyay  (1981)  and  Wei  and  Willmarth  (1989)  investigated  the 
effect  of  Re  on  coherent  structures  in  the  near-wall  region.  Wei  and  Willmarth  (1989)  concluded 
that  as  Re  increases,  the  hairpin  vortices  become  smaller  with  increased  vorticity  in  the  cores.  This 
increased  stretching  of  hairpin  vordees  with  increasing  Re  caused  a  Reynolds  number  effect  on  the 
near-wall  (y+  =  15)  normal  velocity  fluctuations  that  was  not  accounted  for  with  an  inner  variable 
scaling.  As  the  Reynolds  number  (based  upon  the  frec-stream  velocity  and  one-half  channel 
width)  increased  from  3,000  to  40,000,  the  outer  length  scale  8  remained  constant  and  the  inner 
length  scale  v/ux  decreased  by  an  order  of  magnitude.  Such  a  Reynolds  number  effect  could  cause 
the  required  pressure  sensor  diameter  d+  to  decrease  by  an  amount  greater  than  simply  the  decrease 
in  v/ut  •  An  experiment  in  which  Rq  could  be  varied  over  a  wider  range  (1 ,000  <  Re  <  100,000) 
with  pressure  sensors  of  varying  diameter  is  required  to  establish  this  effect.  Current  sensor 
technology  is  inadequate  for  providing  sufficiently  small  sensors  for  such  an  investigation. 


DERIVATION  OF  A  WAVENUMBER  FILTER  DESIGN  FOR 
MEASUREMENTS  AT  HIGH  FREQUENCIES 

In  order  to  design  a  wavenumber  filter,  we  require  a  model  of  the  true  wall  pressure  field. 
Although  we  wish  to  measure  the  two  dimensional  spectrum  d>(ki,0)),  we  must  use  a  model  for  the 
three-dimensional  spectrum  O(ki,k2,co)  in  order  to  account  for  the  two-dimensional  spatial 
response  of  the  sensors.  The  Corcos  (1963)  model  of  the  cross- spectrum  of  the  wall  pressure  field 
is  given  as 

<!>($, T1,<D)  =  a>(co)A(co^Uc)e-‘(“VUc)B(0)ri/Uc) ,  (1) 

where  £,  and  r|  are  the  streamwise  and  spanwise  separation  variables,  respectively.  Experimental 
data  support  exponential  forms  for  A(a)^AJc)  and  B(cor]/Uc)  given  by 

A(co£/Uc)  =  eal“W »  B(cori/Uc)  =  eP  '“W .  (2) 

Farabee  and  Casarella  (1991)  noted  that  the  streamwise  decay  constant  a  reported  from  various 
investigations  varies  from  -0.10  to  -0.19,  and  is  a  function  of  various  flow  parameters.  A  value 
for  a  of -0.125  fits  the  recent  data  of  Farabee  and  Casarella  (1991)  and  also  Keith  and  Barclay 
(1993)  at  the  higher  frequencies.  A  constant  of  -0.125  implies  that  a  90-percent  reduction  in 
A(coE/Uc)  occurs  due  to  the  spatial  decay  of  coherent  structures  as  they  convect  over  three 
wavelengths  in  the  streamwise  direction.  At  lower  values  of  g^/Uc>  similarity  scaling  breaks 

down,  as  shown  by  Farabee  and  Casarella  (1991)  and  Keith  and  Barclay  (1993).  The  convective 
wavenumbers  in  this  range  exhibit  a  rapid  loss  of  spatial  coherence  and  lower  convection 
velocities.  The  similarity  scaling  is  therefore  only  valid  over  a  limited  range  of  convective 
wavenumbers  and  related  frequencies.  Following  Willmarth  and  Roos  (1965),  the  spanwise  decay 
constant  is  taken  as  |3  =  -0.7. 

The  wavenumber-frequency  spectrum  d>(ki,k2,(0o)  is  defined  as 

d>(ki,k2,(Uo)  =  l/4rc2  J  J  <!>({;, Ti dq  .  (3) 

— oo  — oo 

Using  equations  (1)  and  (2),  equation  (3)  was  evaluated  with  the  symbolic  manipulation  program 
MACSYMA  (Rand,  1984)  to  obtain 


3 


<D(ki,k2,Wo)  =  <J>(co0)Fi(ki,kc)F2(k2,kc) , 


(4) 


where 

FKkj.kc)  =  { 1/2  jt  kc)  {-2a  /{a2  +  (kj/kc  +1)2}} 

(5) 

F2(k2,kc)  =  { 1/2 it  kc}  (-2(5  /(p2  +  (k2/kc)2)}  . 

In  transforming  <t>(i;,Ti,coo),  kc  =  ®o/Uc  is  taken  as  constant,  and  the  variation  of  Uc  with  ^  is 
neglected.  The  Corcos  model  of  the  two-dimensional  wavenumber-frequency  spectrum  O(ki,«o) 
is  then 


d>(ki,<o0)  =  J  a>(ki,k2,co0)dk2  = '4>(to0  JF^ki.kc) .  (6) 


The  model  of  the  wavenumber-frequency  spectrum  given  by  equation  (6)  is  shown  in 
dimensional  form  in  figure  1.  The  low  and  high  frequencies  chosen  correspond  to  the  recent  water 
tunnel  investigation  of  Keith  and  Barclay  (1993)  (see  table  1).  There,  the  lowest  measurable 
frequency  was  determined  by  the  acoustic  noise  in  the  facility,  and  the  highest  frequency  by  the 
electronic  noise  from  the  sensor  preamplifiers.  The  cross-spectral  measurements  of  Keith  and 
Barclay  (1993)  agreed  well  with  the  Corcos  (1963)  model  given  by  equations  (1)  and  (2),  with  the 
decay  constant  a  =  -0.125.  The  decrease  of  53  dB  in  the  convective  ridge  peak  level  from  the  low 
to  the  high  frequency  reflects  the  36-dB  decrease  in  the  measured  values  of  <X>(t Oo),  and  the  17-dB 
decrease  in  the  function  Fi(ki,kc).  The  spectrum  O(ki,o>o)  'n  figure  1  may  be  viewed  as  an 
estimate  of  a  typical  spectrum,  based  upon  cross-spectral  measurements  that  could  be  measured 
directly  with  a  wavenumber  filter.  The  high-frequency  portion  of  the  spectrum  will  determine  the 
minimum  sensor  size  and  spacing  requirements  due  to  the  high  wavenumbers  involved.  The 
required  length  of  the  wavenumber  filter  will  increase  significantly  for  the  low-frequency  case  in 
order  to  resolve  the  narrower  convective  ridge. 

# 

We  first  estimate  the  highest  nondimensional  frequency  at  wliich  measurements  of  <J>(ki,coo) 
are  to  be  made  that  are  valid  for  any  turbulent  boundary  layer.  Farabee  arid  Casarella  (1991) 
concluded  that  for  nondimensional  frequencies  cov/ux2  >  0.3,  contributions  to  <J>(co)  are  from 
velocity  fluctuations  in  the  buffer  region,  and  <t>  (co)  scales  with  inner  variables.  We  therefore 
define  ki*=  kjv/ux ,  k2*=  k2v/ux  ,  o>o*=  coov/ux2  ,  Uc*=  Uc/ux ,  k<;*=  coq */Uc*=  kev/ux  , 
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Figure  1 .  Dimensional  Wavenumber-Frequency  Spectra  at  High  and  Low  Frequencies  Estimated 
From  the  Cross-Spectral  Measurements  of  Keith  and  Barclay  (1993) 


Table  1.  Dimensional  Wavenumber  Filter  Parameters  for  Measurements  at  the 
High  Frequencies  in  Selected  Turbulent  Boundary  Layers 
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figure  3.  We  have  used  rectangular  rather  than  circular  sensors  to  simplify  the  analysis.  The 
sensor  voltage  outputs  p(nAx*,t*),  n  =  0,  1,2, N-l,  are  acquired  simultaneously  in  time.  The 
analog  data  are  first  bandpass  filtered  in  the  time  domain.  This  may  be  viewed  as  a  two-step 
process.  Unwanted  acoustic  energy  inherent  to  the  flow  facility  (typically  at  low  frequencies)  is 
removed  with  frequency  filtering.  These  data  are  then  again  filtered  at  the  Nyquist  frequency  to 
avoid  temporal  aliasing.  The  sensor  outputs  are  next  digitized  in  the  time  domain  to  generate 
p(nAx  ,mAt  ).  The  discrete  temporal  Fourier  transform  of  each  sensor  output  is  then 

M-t 

PrfnAx^Wo*)  =  At*/2n  £  wt(m)p(nAx*,rnAt*)e'i(0o*mAl*  .  (10) 

m=0 

Here,  W[(m)  is  the  temporal  window,  and  r  refers  to  the  temporal  record  for  the  ensemble 
averaging,  r  =  1,2, ....  R.  Contemporary  analog-to-digital  converters  and  mass  storage 
capabilities  allow  sufficient  sampling.  At*,  and  record  length,  M,  to  minimize  errors  in 
Pr(nAx*,oa0*).  The  discrete  spatial  Fourier  transform  of  Pr(nAx’t’,0)o*)  is  then 

N-l 

Hr(ki*,co0*)  =  Ax*/2jc  X  wS(n)Pr(nAx*,0)o*)e'iki*nAx*  .  (11) 

n=0 

The  measured  wavenumber-frequency  spectrum  is  thus 

Sr(k1*,G)o*)  =  [l/(NAx*MAt*CiC2)]IHr(ki*,w0*)l2  ,  (12) 

where  Q  and  C2  are  corrections  for  the  temporal  and  spatial  windows  given  by 

M-l 

Cl  =  (1/M)  £  wt2(j), 
j=o 

(13) 

N-l 

C2  =  (1/N)  X  ws2(j)  . 

j=0 

We  choose  the  Hanning  window  given  by  w(q)  =  1  -  cos2[7tq/(Q  -  1)],  q  =  0,  1,  2, ...,  Q-l,  with 
side  lobes  32  dB  below  the  main  lobe,  for  both  the  temporal  and  spatial  cases.  The  ensemble 
averaged  wavenumber-frequency  spectrum  is  then 
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(14) 


R 

S(ki*,(Oo*)  =  1/R  X  Sr(ki*,o)o*l  . 

t-i 

Equations  (10)  through  (14)  may  be  used  directly  to  determine  S(k)*.coo*)-  However,  existing 
fast  Fourier  transform  (FFT)  routines  arc  preferable  so  as  to  compute  the  discrete  Fourier 
transforms  more  efficiently.  A  block  diagram  showing  the  sequence  of  steps  in  the  computation  of 
S(ki  ,C0o  )  is  shown  in  figure  4. 

We  now  consider  the  errors  that  result  from  the  spatial  response  of  the  sensors  and  also  from 
the  wavenumber  filter.  We  note  that  filtering  in  the  spatial  domain  is  analogous  to  filtering  in  the 
time  domain.  However,  aliasing  in  the  spatial  domain  is  more  difficult  to  eliminate,  as  will  be 
shown.  The  relationship  between  the  measured  and  the  true  spectrum  is  given  by 

oo  oo 

S(k!*,C0o*)=  J  J  E(Ki*,K2*,0)o*)1Q*(Ki*,K2*)12 

—AO  —  CO 

(l/NAx*27tC2)IZ*(K1*  -  k1*)l2dK1*dK2*,  (15) 

where  the  normalized  sensor  response  is  given  by 

Q*(Ki*,k2*)  =  [sin(K  1*L1*/2)/(K1*L1*/2)][sin(K2*L2*/2)/(K2*L2*/2)] ,  (16) 

and  the  wavenumber  filter  response  is 
N-l 

ZW)  =  Ax*  £  ws(q)e-iKiNAx*  .  (17) 

q=0 


We  first  consider  errors  resulting  from  the  sensor  response  Q*(Ki’k,K2¥).  The  Corcos  model 
for  E*(k1*,k2*,co0*)  given  by  equation  (7)  and  the  sensor  response  Q*(Ki*,k2*)  given  by 
equation  (16)  are  both  separable  in  Ki*  and  K2*.  We  define 

E*(K1*,K2*,a)o*)lQ*(K1*,K2*)|2  =  g,*(K1*,L1*)g2*(K2*,L2*)  ,  (18) 

where 
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Figure  4.  Signal  Processing  Flow  Chart 


g1*(K1*,L1*)  =  F*,(K1*)[sin(K1*L1*/2)/(K1*L1*/2)l2  , 
S2*(k2*,L2*)  =  F*2(K2*)(sin(K2*L2*/2)/(K2*L2*/2)]2  , 


(19) 


and  F*i(Ki*)  and  F*2(k2*)  are  given  by  equation  (8).  From  equation  (15),  the  effect  of  the 
spanwise  dimension  L2*  on  S(k]*,<Do*)  is  given  by 


G2*(L2*)  =  J  g*2(K2*,L2*)dK2*  .  (20) 

For  a  particular  L2*.  the  measured  spectrum  S(ki*,coo*)  is  reduced  by  a  constant  amount  at  all  ki* 
and  coo*  Equation  (20)  was  evaluated  using  MACSYMA  to  obtain 

G2*(L2*)  =  2[e(P  -  pkc*L2*  -  l]/(|)kc*L2*)2 .  (21) 

The  function  G2*(L2*)  is  shown  in  figure  5.  For  L2*=  2,  a  10-percent  (0.5  dB)  reduction  in  the 
level  of  S(ki*,C0o*)  occurs,  and  for  L2*=  5,  a  20-perccnt  (1  dB)  reduction  occurs.  Further  insight 
into  the  spatial  response  of  the  sensor  may  be  gained  by  examining  the  function  g2*(K2*>E2*). 
This  function  is  shown  in  figure  6  for  L2**  0  (true  spectrum),  5, 10,  20,  and  100.  Here,  both  the 
main  lobe  of  the  spanwise  sensor  response  and  the  dominant  energy  in  the  true  spectrum  are 
centered  around  K  *  -  0,  The  spatial  filtering  of  the  sensor  is  significant  at  higher  positive  and 
negative  wavenumbers  due  to  the  broad  energy  content  of  the  true  spectrum.  The  wavenumber 
range  - 1  <  k2*  <  1  of  figure  6  contains  approximately  90  percent  of  the  energy  existing  in  the  true 
spectrum  over  all  wavenumbers.  For  values  of  L2*  greater  than  5,  the  sensor  filtering  dominates, 
as  evident  from  the  side  lobe  patterns  in  the  figure.  We  choose  L2*=  5  as  an  acceptable  design 
parameter,  but  note  that  smaller  sensors  are  desirable. 

An  acceptable  size  for  Lj*  can  only  be  determined  by  considering  the  sensor  response  in 
conjunction  with  the  wavenumber  filter  response,  as  evident  from  equation  (15).  The  filtering  of 
high  wavenumbers  by  the  sensor  limits  the  energy  that  is  aliased  into  the  wavenumber  range  of 
interest.  We  first  compare  the  effect  of  the  sensor  response  alone  on  the  true  spectrum,  as  this 
determines  the  field  that  the  wavenumber  filter  operates  on.  Since  we  are  measuring  S(ki*,CDo  ), 
the  function  gi*(Ki*,Li*)  itself  determines  the  filtering  effect  of  the  sensor.  This  function  is 
shown  in  figure  7,  for  Li  =0  (true  spectrum),  5,  10,  20,  and  100.  The  main  lobe  of  the  sensor 
response  is  again  centered  around  Ki*=  0,  but  the  dominant  energy  in  the  tr  ue  spectrum  is  centered 


12 


G'i(Li').  Eqn  (a4) 
G'2(l2').Eqn  (21) 
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Figure  6.  The  Product  of  the  Spanwise  Sensor  Response  (for  0  <.  L2*  <!  100)  and 
the  True  Spectrum  for  the  High-Frequency  Case 


Figure  7.  The  Product  of  the  Streamwise  Sensor  Response  (for  0  £  Li*  £  100)  and 
the  True  Spectrum  for  the  High-Frequency  Case 

around  Ki*  =  -0.2.  For  Li*=  5.  the  peak  in  gi*(Ki*,Li*)  occurs  at  Ki*=  -0.2  ,  and  the  peak  level 
is  reduced  by  8  percent  (or  0.4  dB).  For  Li*  =  10,  the  peak  is  shifted  by  0.5  percent  toward  the 
origin,  and  the  peak  level  is  reduced  by  29  percent  (or  1.5  dB).  For  larger  values  of  Lj*.  the 
sensor  filtering  dominates  the  measured  spectrum,  as  shown  by  the  side  lobe  patterns.  We  select 
Lj*  =  5,  and  determine  the  effect  of  the  wavenumber  filter  on  the  measured  spectrum. 

With  L-2*  =  5,  a  20-percent  reduction  was  shown  to  occur  in  S(ki*,coo*).  We  express 
equation  (15)  as 


S(ki*,Oi0*)  =  C3  J  g,Vi*)(l/NAx*2jc  C2)iZ*(iCi‘  - 


m\ 

\^9 


where  C3  =  0.8.  In  order  to  determine  die  effect  of  the  wavenumber  filter  on  the  measured 
spectrum,  we  consider  the  quantity  S(ki  ,u>o  )/C3  defined  as  Si  (ki  ,C0q  )  and  do  not  include  the 
20-percent  reduction,  which  is  constant  at  all  ki*.  The  filter  function  given  by 


V(Ki*  -  kj*)  =  (l/NAx*2jtC2)iZ*(Ki*  -  ki*)P 


(23) 
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is  periodic  in  wavenumber  at  period  2tc/Ax  .  The  number  of  sensors  N  and  spacing  Ax 
determine  the  wavenumber  bin  width  Ak*  =  2rc/NAx*.  We  initially  define  Ax*  =  10  (allowing  for 
finite  spacing  between  the  sensors)  and  N  =  100.  The  filter  function  V(K'i*  -  k|*)  is  shown  in 
figure  8,  where  the  main  lobe  is  centered  at  ki  =  -0.2.  In  comparison,  we  also  show  the  filter 
function  for  the  case  Ax*  =  5  and  N  =  100  with  the  main  lobe  centered  at  ki  =  -0.6  (for  clarity). 
Decreasing  the  spacing  from  10  to  5  viscous  lengths  (with  the  number  of  sensors  held  constant) 
increases  the  wavenumbers  at  which  the  aliasing  lobes  occur  by  a  factor  of  two.  In  addition,  the 
main  lobe  (and  aliasing  lobes)  becomes  slightly  broader,  and  the  wavenumber  bin  width  increases 
by  a  factor  of  2.  It  is  therefore  desirable  to  maximize  the  number  of  sensors  to  achieve  a  narrow 
main  lobe  and  bin  width.  In  addition,  minimizing  the  sensor  spacing  causes  the  aliasing  lobes  to 
occur  at  wavenumbers  where  the  spectral  levels  are  low.  We  also  utilize  the  filtering  provided  by 
the  sensors  (previously  discussed)  to  reduce  the  spectral  levels  at  which  the  aliasing  lobes  occur, 
since  the  field  that  the  filter  function  operates  on  is  effectively  prefiltered  by  the  sensors. 

The  measured  spectrum  Si(ki*,coo*)  was  computed  numerically  for  three  different  cases,  as 
shown  in  figure  9.  The  first  case  is  for  100  infinitesimally  small  sensors  (Li*  =  0)  and  Ax*  =10  in 
order  to  illustrate  the  effect  of  the  filter  function  V(K]*  -  kt*)  alone  on  the  measured  spectrum. 

The  measured  spectrum  is  periodic  at  2rc/10,  and  agrees  well  with  the  true  spectrum  in  the 
convective  ridge  region  for  -0.3 1  <  ki*  <  -0.09.  For  kj*  >  -0.09,  the  measured  levels  are  higher 
than  the  true  levels  due  to  aliasing.  The  second  case  is  for  100  sensors  with  Lj  =5  and  Ax  =10. 
The  measured  levels  are  slightly  lower  than  the  true  levels  for  ki*  <  -0.2  due  to  the  filtering  from 
the  sensors  (note  the  measured  levels  for  case  2  and  case  3  are  identical  over  this  wavenumber 
range).  For  -0.2  <  ki*  <  0,  the  measured  levels  agree  well  with  the  true  levels,  with  the  aliased 
energy  reduced  due  to  the  sensor  filtering.  For0<ki  ,  aliased  energy  is  apparent.  The  third  case 
involves  200  sensors  with  Li  =5  and  Ax  =  5.  The  measured  levels  agree  well  with  the  second 
case  for  ki  <  0.  For  0  <  ki  the  measured  levels  agree  very  well  with  the  true  spectrum.  The 
aliased  energy1  apparent  for  cases  1  and  2  over  this  wavenumber  range  is  reduced.  This  is  due  to 
the  smaller  sensor  spacing,  which  places  the  aliasing  lobes  farther  from  the  main  lobe.  Here,  the 
smaller  spacing  leads  to  increased  accuracy  at  low  values  of  lki*l.  We  select  the  parameters  for  the 
second  case,  which  result  in  a  good  estimate  of  the  true  spectrum  around  the  convective  ridge  and 
allow  for  finite  spacing  between  the  sensors,  which  is  desirable  for  fabricating  the  array  of 
sensors. 
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Figure  8.  The  Wavenumber  Filter  Responses  Compared  With  the  True  Spectrum 

for  the  High-Frequency  Case 


DATA  POINTS  ARE  DISPLAYED  AT  EVERY  SECOND  WAVENUMBER  BIN  FOR  CLARITY 


Figure  9.  Comparison  of  the  Measured  Spectrum  (for  Three  Different  Wavenumber 
Filters)  With  the  True  Spectrum  at  the.  High  Frequency 
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MEASUREMENTS  AT  LOW  FREQUENCIES 


We  now  consider  the  measurement  of  the  wavenumber  frequency  spectrum  at  low  frequencies. 
The  lowest  frequency  at  which  measurements  may  be  made  is  often  determined  by  the  background 
acoustic  noise  in  the  experimental  facility.  The  use  of  noise  cancellation  techniques  as  discussed 
by  Farabee  and  Casarella  (1991)  can  allow  measurements  to  be  made  at  frequencies  where  acoustic 
noise  exists.  However,  an  absolute  low-frequency  limit  valid  for  any  boundary  layer  cannot  be 
specified.  With  regard  to  the  turbulence  physics,  it  is  desirable  to  make  measurements  at  as  low  a 
frequency  as  is  physically  possible  for  a  particular  facility  or  test.  Here  we  arbitrarily  choose  a  low 
frequency  and  determine  suitable  parameters  for  the  measurement  of  the  wavenumber  frequency 
spectrum. 

Farabee  and  Casarella  (1991)  defined  a  frequency  region  5  <,  coS/ux  £  100,  where  the 
dominant  sources  for  O(w)  are  from  velocity  fluctuations  in  the  outer  region  of  the  turbulent 
boundary  layer.  We  choose  (0o8/ux  =  20  which  falls  within  this  range.  The  relationship  between 
the  outer  and  inner  frequency  scaling  is  given  by  o>o*Rx  =  to05/ux  .  where  Rx  =  5ux/v  .  With  Rx 
=  2010  (for  the  highest  Reynolds  number  of  Farabee  and  Casarella  (1991)),  too*  =  .01.  The  value 
for  the  streamwise  decay  constant  a  was  taken  as  -0.24,  from  the  results  of  Farabee  and  Casarella 
(1991)  at  the  lower  frequencies.  The  spanwise  decay  constant  P  was  taken  as  -0.7,  from  the 
results  of  Willmarth  and  Roos  (1965).  For  convenience,  we  will  retain  the  inner  variable  scaling 
for  displaying  the  spectra,  but  note  an  outer  variable  scaling  is  appropriate  at  low  frequencies,  as 
discussed  by  Farabee  and  Casarella  (1991)  and  Keith  et  al.  (1992). 

For  measurements  at  low  frequency,  sensor  size  requirements  are  significantly  less  restrictive 
than  for  those  at  high  frequency.  This  is  due  to  the  more  narrow  distribution  of  energy  in  both  ki* 
and  k2  in  the  true  spectrum  at  low  frequencies,  as  shown  in  figures  2  and  10.  Larger  values  of 
Li*  are  in  fact  desirable  for  reducing  aliased  energy,  as  previously  discussed.  The  reduction  in 
spectral  levels  due  to  the  spanwise  length  is  negligible  (within  2  percent)  for  L2*  ^  100,  as  shown 
by  the  values  for  G2*(L2*)  in  figure  11.  We  define  L2*  =  100,  but  note  that  larger  values  are 
acceptable. 

A  smaller  wavenumber  bin  width  is  required  to  resolve  the  narrower  distribution  of  energy 
around  the  convective  ridge.  The  overall  length  of  the  wavenumber  filter  must  therefore  be 
increased.  We  define  the  bandwidth  Ak  over  which  the  spectrum  levels  are  within  3  dB  of  the 
peak  level  al  the  convective  ridge.  For  the  high-frequency  case,  Ak  contains  eight  wavenumber 
bins.  We  choose  100  sensors  with  length  Lj*  =  100  and  spacing  Ax*  =  1480,  which  yields  eight 
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Attenuation 


Figure  10.  Nondimensional  Wavenumber-Frequency  Spectrum  for  the  High-  and 
Low-Frequency  Cases  From  the  Corcos  Model 
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Figure  1 1.  Spectral  Attenuation  at  die  Low  Frequency  Due  to  the  Spanwise  or 

Streamwise  Sensor  Length 
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Figure  12.  Comparison  of  the  Measured  Spectrum  (for  Three  Different  Wavenumber 
Filters)  With  the  True  Spectrum  at  the  Low  Frequency 

bins  over  Ak  for  the  low-frequency  case.  The  measured  spectrum  for  this  case  is  shown  in  figure 
12.  The  convective  ridge  region  is  accurately  resolved,  and  the  increased  levels  at  higher  and 
lower  wavenumbers  occur  from  aliasing.  The  overall  length  of  the  wavenumber  filter  for  this  case 
is  148,000  viscous  lengths.  Also  shown  in  figure  12  are  the  cases  for  Li*  =  500  and  Li*  =  1000, 
with  the  same  number  of  sensors  and  sensor  spacing.  Increasing  the  sensor  length  increases  the 
wavenumber  range  over  which  the  convective  ridge  region  is  accurately  resolved,  due  to  the  filter- 
ing  of  wavenumbers  that  are  aliased.  For  the  Li  =  1000  case,  the  spectral  levels  are  below  the 
true  levels  at  the  higher  negative  wavenumbers,  but  agree  well  at  the  higher  positive  wavenumbers. 
The  Li*  =  500  length  is  chosen  for  the  purpose  of  designing  a  filter  at  this  low  frequency. 

Clearly,  separate  wavenumber  filters  are  required  at  high  and  low  frequencies.  To  achieve  an 
overall  filter  length  of  148,000,  required  for  tire  low-frequency  case  (as  shown),  the  number  of 
sensors  needed  using  a  10  viscous  length  spacing  (corresponding  to  the  high-frequency  filter)  is 
14,800.  This  is  impractical  in  view  of  the  signal  processing  requirements  as  well  as  the 
wavenumber  filter  fabrication.  For  a  particular  application,  we  may  need  to  design  a  number  of 
filters  to  cover  the  frequency  ranges  of  interest.  The  dynamic  range  required  for  individual  sensors 
therefore  depends  on  the  frequency  range  of  the  wavenumber  filter  for  which  they  will  be  used. 
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DIMENSIONAL  YVAVKN1IMHKU  MLTLU  I'AUAMI  ll  US  Kilt 
TUKHULLNT  HOUNDAKV  LAYl'.ltS  IN  AIK  AND  \YAi»  U 


Wc  consider  boimduty  layets  from  paitleulat  wind  (imm.'l  mul  wuict  uinnvl  invcMigniionv  n<*i1 
also  an  external  boundary  layci  that  could  exist  tit  tut  mulct  son  upplU  ubmi,  mu  h  u\  How  ov.  ■  « 
hull.  The  boundary  layci  paiuntelcis  and  the  pto|toscd  wavcmtmU  i  Itliu  piitunu tvis lm  tin  Ingli 
frequency  ease  arc  given  in  table  1.  Keguiding  the  two  wind  tunnel  iitvcMlpitiit »n%.  Wtllmnitlt  m,.t 
Wooldridge  (19(>2)  worked  at  a  high  ltcc-Micum  velocity,  (ml  tlteli  till,  ket  Utuiiilniv  Into  l>  <1 1"  a 
viscous  length  scale  comparable  to  that  ot  hatabce  and  t'usuirllu  t  I'^M  ).  1  ot  tin  tn\<  \iig,u...n- 
using  wuler,  the  boundary  luyet  thicknesses  wete  enmptuuhle.  hut  the  hlghei  I  tee  Mourn  v.'l<*  nv 
of  Carey  ot  al,  (,1907)  led  to  a  sntuller  viscous  length  sealc  than  that  <>l  Keith  mid  linn  l«\  i  I’hmi 
l;or  the  undersea  application,  the  Ixmnduiy  laye.t  putuniclm  wete  uhtained  f i«*m  t  ’»«!•  1 1 •  > *.  i. 

ideal  sniootli  wall,  zero  pressure  gradient,  tlat  plate  case.  The  likudon  limn  tin  tmlml.  ni  U.im.i  ns 
layer  otigitl  wus  taken  us  uppioxitnutely  UK)  It.  The  viscous  length  m  air  is  i  mtipmuhl*'  to  tin  mlm 
water  ilows,  and  the  boundary  layer  is  thicket. 

The  required  sensor  sizes  ior  these  eases  ate  on  the  mdei  ot  It)  miemnv  1  In  n  qulnd  l.-ngd. 

NAx  of  die  wuvatumber  iilter  (based  upon  1(X)  seusois  at  ll)  viscous  length  spin  tugt  \no>  >  . . . 

the  order  of  0.1  to  1.0  cm.  Thcd*  vahtesuctuallyiisedtucaihlnvestigutionatvrtUogK.il  tin 
small  viscous  length  scules  for  the  wutei  ilows  led  to  huge  d*  values.  I  ot  the  wain  ninm  I 
investigation  of  Keith  and  lhuclay,  tlie  actual  setisin  illumetei  was  0  ?  cm,  will.  It  i  mn  spomt.  .1 1. . 
444  viscous  lengths.  The  requited  sensor  sizes,  us  well  as  the  seusm  upm  lugs,  tm  tln  o  .  mo 
emphasize  tin:  lequiiement  ior  mietoialnicaiioit  technology  lm  such  uppln  atimts,  pmtn  ulmh  m 
water. 

CONCLUSIONS  AND  HIst’OMMLNDATlONS 

Hie  analysis  presented  here  shows  that  the  tneiisuietueni  oi  tin  Indent  wall  piesstue 
wuvenumlxM'-frcquency  speeua  at  high  ftequeneles  lot  Ky  values  ol  ptui  deal  tuteiest  m  ,m  m  w.it.‘i 
requires  microfubi ieation  technology.  The  spauwisc,  as  well  as  the  stieamwise,  seusm  length  is 
significant  in  reducing  the  measured  spccual  levels.  I'm  the  measuiemeut  ol  the  wuwiunnh.  i 
frequency  siieenum  <I>(ki,<o),  the  stieamwise  seusm  length  also  limits  the  s.-nsm  sp.n  tug  mnl 
affects  the  amount  of  aliased  energy, 

We  have  considered  two  particular  frequencies  tei]uuiug  two  waveiiumbi  i  tih<-i  designs  At 
high  frequencies,  a  short  overall  filter  length  and  small  seusm  size  and  sntsm  •  p.n  mg  m.  i.  <|Uii.  d 
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\*  low  hv>|n«'m  iek.  ilir*  hlivi  length  iMvvn  nuli’is  ol  iiiii|*nitiulc*  greater,  ns  is  the  sensor  spacing. 

In  •ililMimi,  m’iikiiv  (luu  me  iwii  oi«loi  v  ill  magnitude  potter  in  stieaimvi.se  length  aio  required  to 

f«*lll|,  •  «ll«li>il  1'tll‘l  J.V 

\N »>  lmv«>  imumeit  ideal  it'i  tiuigului  sensois  with  known  wiivciiiiuiIki  frequency  responses.  In 
i«i*»  in  ♦».  i'wh  ilnmgh  iln  luijiit'in  y  lrsponse  ol  sensois  can  be  nicnsuicd,  (lie  wuvcuumbcr- 
iivi|iieui  v  n  *|>.«iiii'  iii  unknown  ami  noi  easily  mcasiued.  While  icctuugular  scnsois  may  l>c 
pi»leMilile  ini  Miiunilm  lining  w<'l<  mil  iulnluh  uilmi  technology,  out  unulysls  can  also  l>c  extended 
in  «  iii  mUi  winkoM  We  have  uvnl  the  Tokos  model  line,  hut  olhet  models  may  Ik  considered. 
NhiIh  i  tin-  mu  eitaimh  s  tine  to  the  settsoi  u‘ spouses  not  the  selected  model  will  slgililicuiltiy  uliCCt 
tlw>  i  mit  liiMom  pu  kiMileil  tit  lilts  pupei . 

Ai  ill*-  pu  »*ni  lime,  the  an  man  oieasuiement  ol  ensemble  uvomgal ' vavenumber Ircqncticy 
«pi*»  n a  xi  vipillilmiim  Uuntihny  luyei  wall  piesstue  lluctuallouii  is  ol  sigull  leant  interest  because  ol 
lit#  U>  k  ol  |Hili|ulu«l  ihiiit  I  lowevei,  the  liiteipietuiioii  ol  such  measurements  in  rclution  to  the 
pits  in  ml  ih#  tin  hull  ut  hiiumhuy  layei  will  icqulio  eomliiloual  sampling;  a;  well  as  parametric 
iiiithi'k  I  he  i-lli  1 1  ol  Hu  on  the  measuivd  s|Keliunt  will  Ik  ol  piuticular  interest,  especially  lor 
i  iiinpmikuii  mill  iv kill o  imm  uiimmcal  slmule.lloii.s  that  me  coiistiuined  to  lower  Ro  values,  hi 
•nil  Inn  mi,  the  iiuitMueil  k|Ki  ini  imm  |kiui>IkiI  01  uoiicquililuiuin  Uiuudaiy  layers  will  also  provide 
ridilltioiiril  phy  kii  nl  insight  It  Is  ho|Vil  thill  mlctohibi  icaiioii  technology  will  allow  such 
mi  iiktiii  limits  to  Ik  pm  nil  d  In  the  Inline 
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APPENDIX 

SENSOR  SIZE  EFFECTS  ON  MEASURED  FREQUENCY  SPECTRA 

Although  we  have  focused  on  the  measurement  of  the  wavenumber-frequency  spectrum, 
microfabrication  technology  is  important  for  single-point  measurements  of  the  frequency  spectrum 
Oml Wo)-  For  our  case,  the  relationship  between  the  measured  and  true  (one  dimensional) 

frequency  spectrum  may  be  expressed  as 

oo  oo 

Om(C0o)=<I>l(0)o)  j  J  gl*(K1*,L1*)g2,t‘(K2*,L2*)dK1*dK2*  ,  (A-l) 

_ oo  .00 

which  reduces  to 

Om(o>o)  =  ,  (A-2) 

where  G2*(L2*)  is  given  by  equation  (21)  and 


Gt*(L!*)=  J  g*1(K1*,L1*)dK1^ 


Equation  (A-3)  was  evaluated  using  MACSYMA  to  obtain 

G!*(Li*)  =  2[e(a  kc*Ll*)(2asin(kc*L1*)  +  (a  2  -  l)cos(kc*L1*)) 

-  a(a2  +  l)kc*Li*-  a2  +  l]/f(a  4  +  2a  2  +  l)kc*2L!*2] .  (A-‘ 

The  functions  Gi*(Li*)  and  G2*(L2*)  are  shown  in  figure  5  for  the  high-frequency  case.  From 
equations  (21)  and  (A-4), 


dG2*/dL2*|  L2*=q  =  5  dGi^/dLi*  L*sQ. 


For  small  values  of  Lj  and  L2  ,  the  spanwise  sensor  length  therefore  causes  a  greater  reduction  in 
spectral  level  than  the  same  streamwise  length.  For  Lj*  =  L2*  =  10.4,  the  spanwise  and 
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streamwise  effects  are  comparable.  For  larger  values  of  Lj*  ard  L2*,  tire  streamwise  length 
dominates  the  reduction. 

The  functions  Gi*(Li*)  and  G2*(L2*)  are  shown  in  figure  1 1  for  the  low-frequency  case. 
These  functions  are  similar  to  the  high-frequency  case,  with  the  L2*  effect  dominant  for  sensor 
lengths  of  less  than  2500,  and  the  Li*  effect  dominant  at  larger  values.  However,  the  actual  values 
of  Li*  and  L2*  required  for  a  particular  reduction  are  generally  two  orders  of  magnitude  greater 
than  for  the  high-frequency  case.  This  reflects  the  change  in  the  energy  distribution  in  the  true 
spectrum  at  low  frequencies,  as  shown  in  figures  2  and  10.  The  reductions  are  negligible  for  Li* 
and  L2  <,  100.  For  measurements  at  low  frequency,  sensor  size  requirements  are  therefore 
significantly  less  restrictive  than  for  those  at  high  frequency. 


A-2 


INITIAL  DISTRIBUTION  LIST 


Addressee 

Defense  Technical  Information  Center 

Northwestern  University,  Dept,  of  Mechanical  Engineering  (Prof.  R.  Leuptow) 

University  of  Connecticut,  Dept,  of  Mechanical  Engineering 
(Prof.  J.  Bennett,  Prof.  L,  Langston) 

OASIS  (Dr.  C.  Gedney) 

Rutgers  University,  Dept,  of  Mechanical  and  Aerospace  Engineering 
(Dr.  T.  Wei) 

U.S.  Coast  Guard  Academy  (Dr.  V.  Wilczynski,  K.  Colella) 

University  of  Minnesota,  Twin  Cities  Campus,  Dept,  of  Aerospace  Engineering 
and  Mechanics  (Dr.  A.  Alving) 

Univ  'tsity  of  Minnesota,  Twin  Cities  Campus,  Dept,  of  Mechanical  Engineering 
(Prof.  T.  W.  Simon) 

Office  of  Naval  Research  (Dr.  L.  P.  Purtell,  Dr.  E.  P.  Rood,  Dr.  K.  Ng) 
NASALRC  (Dr.  J.  C.  Hardin) 

The  Catholic  University  of  America  (Prof.  M.  Casarella) 

Tufts  University  (Prof.  V.  P.  Manno) 

Naval  Surface  Warfare  Center,  Carderock  Division  (Dr.  T.  Farabee,  P.  Chang) 
University  of  Michigan  (Dr.  P.  B.  Washabaugh) 

TSI,  Inc.  (G.  Reynolds) 

General  Atomics  (Dr.  C.  R.  Harder,  P.  W.  Trester) 


No,  of  Copies 

12 

1 

2 

1 

1 

2 

1 

1 

3 

1 

1 

1 

2 

1 

1 

2 


